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ABSTRACT 

We report the analysis of the rotational properties of our sample of Herbig Ae/Be (HAeBe) 
and related stars for which we have obtained high-resolution spectropolarimetric observa- 
tions. Using the projected rotational velocities measured at the surface of the stars, we have 
calculated the angular momentum of the sample and plotted it as a function of age. We have 
then compared the angular momentum and the v sin i distributions of the magnetic to the non- 
magnetic HAeBe stars. Finally we have predicted the vsin/ of the non-magnetic, non-binary 
("normal") stars in our sample when they reach the ZAMS, and compared them to various 
catalogues of the v sin i of main-sequence stars. First, we observe that magnetic HAeBe stars 
are much slower rotators than normal stars, indicating that they have been more efficiently 
braked than the normal stars. In fact, the magnetic stars have already lost most of their an- 
gular momentum, despite their young ages (lower than 1 Myr for some of them). Secondly, 
our analysis suggests that the low mass (1.5 < M < 5 M ) normal HAeBe stars evolve with 
constant angular momentum towards the ZAMS, while the high-mass normal HAeBe stars 
(M > 5 M ) are losing angular momentum. We propose that winds, which are expected to be 
stronger in massive stars, are at the origin of this phenomenon. 
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1 INTRODUCTION 

Among intermediate mass A/B stars on the main sequence 
(MS), a sub-group called chemically peculiar Ap/Bp stars shows 
very specific properties : abundance anomalies at their surface, slow 
rotation (rotation periods longer than 1 day), and strong magnetic 
fields. While the abundance anomalies are believed to result from 
diffusion within their surface layers due to the com petition betwee n 
radiative levitation and gravitational settling (e.g. Micha udl l 19700 . 
the origin of the slow rotation is less well understood. The me- 
chanism inducing the slow rotation of Ap/Bp stars is likely related 
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to their strong magnetic fields. Stepien| i feOOOh discussed different 
theories, and he concludes that magnetic braking must occur du- 
ring the pre-main sequence (PMS) phase in order to reproduce the 
angular momentum observations of the MS Ap/Bp stars. Stepien 
demonstrated that magnetic coupling of a PMS star with its accre- 
tion disk would slow the star and reduce its rotation period to a few 
days. In order to produce the very slowest rotators, with observed 
rotation periods greater than a month, the disk must disappear suf- 
ficiently early during the PMS phase, allowing strong magnetised 
winds to carry away a large quantity of angular momentum until 
the stars reach the zero-age main-sequence (ZAMS). 

The PMS progenitors of MS A/B stars are found among Her- 
big Ae/Be (HAeBe) stars. These stars have emission lines in their 
spectra, lie in an obscured region, illuminat e fairly bright nebulae, 
and show a pronounce d infrared excess |Herbig 196Q; iThe et all 
1 1994 IVieiraetal .1120031) . all of which are the result of their young 
age. All of them do not necessarily show all Herbig characteris- 
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tics, but all of them have infrared excess with an abnormal extinc- 
tion law (compared to classical Be stars). The PMS st atus of these 
stars w as first confirmed by the spectroscopic study of lStrom et al.l 
(1972), by showing that their surface gravities are systematically 
lower than those of their MS counterparts. 

Until recently we had very few observational constraints on 
the magnetic fields and rotation of these objects. To our know- 
ledge two thorough observational studies of the evolution of an- 
gular momentum o f intermediate-mass sta rs during the PMS phase 
have been realised. lB6hm & Catafddl995l) . using MUSICOS high- 
resolution spectra, performed a statistical study of the projected ro- 
tational velocities (v sin f) of 29 HAeBe stars. They studied the evo- 
lution of the spin angular momentum of intermediate mass stars du- 
ring the PMS phase by comparing the v sin i of field HAeBe stars to 
that of MS A/B stars in young open clusters. They concluded that, 
assuming solid body rotation within the star, the evolution depends 
on the mass of the star : low mass HAeBe stars (M < 2.6 M Q ) seem 
to lose a large fraction of their angular momentum during their 
PMS evolution ; intermediate mass HAeBe stars (2.6 < M ^ 4 M Q ) 
seem not to lose a significant fraction of their angular momentum ; 
and higher mass HAeBe stars (M ^ 4 M Q ) seem to gain a large 
amount of angular momentum. However, if the stellar internal ro- 
tation varies as r~ 2 , the observations of HAeBe stars and A/B stars 
in young clusters are consistent with conservation of total angular 
mo mentum at al l masse s. 

IWolff et~aIU2004h studied the v sin i of 145 stars between 0.4 
and 14 M Q located in the Orion star-forming complex. They argued 
that observations of the surface rotation velocities at different evo- 
lutionary stages of the PMS phase can provide information about 
the internal angular momentum profile of stars. They therefore dis- 
tinguished between stars still on their convective track (the Hayashi 
phase) and more evolved stars on their radiative track. Their conclu- 
sions lead to the following scenario concerning intermediate mass 
stars (M > 1.5 M Q ) : stars lose angular momentum before they start 
the PMS phase, a decoupling of the angular momentum seen at the 
surface from the angular momentum in the interior occurs when the 
stars make the transition from convective to radiative PMS evolu- 
tion, after which the angular momentum is conserved during the 
PMS radiative evolution. 

Both analyses gave very interesting results, and should now 
be discussed in the framework of a scenario of angular momen- 
tum evolution that includes magnetic fields. With this aim, we have 
used the large survey of 70 HAeBe stars that we performed using 
the high-resolution spectropolarimeters ESPaDOnS at the Canada- 
France-Hawaii Telescope (CFHT, Hawaii), and Narval at the Teles- 
cope Bernard Lyot (T BL, Pic du Midi, France). In a previous paper 
( I Alecian et ai]|2012i . paper I), we have described the whole obser- 
ved sample, including the measurement of their projected rotation 
velocities (vsin/). We now present the second paper of the series 
describing the analyses of the rotational velocities and angular mo- 
mentum evolution of the sample. Section 2 gives a short summary 
of the observations and the measurements reported in paper I. Sec- 
tion 3 reports our findings concerning the angular momentum evo- 
lution among our sample, while Section 4 describes the analysis 
performed on the projected rotational velocities. Finally, in Section 
5 we conclude and discuss our results. 



2 SUMMARY OF THE SAMPLE AND OBSERVATIONS 

We used the ESPaDOnS and Narval instruments in polarime- 
tric mode to obtain 132 spectra of 70 HAeBe stars with a resolving 
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Figure 1. Magnetic (red squares) and non-magnetic (black points) Herbig 
Ae/Be stars plotted in a Hertzsprung-Russell diagram. Open circles corres- 
pond to HD 98922 (above the birthline) and IL Cep (below the ZAMS) 
that fall outside of the PMS region of the HR diagram, whose positions 
cannot be reproduced with the theoretical evolutionary tracks considered 
in this paper. The CESAM PMS evolutionary tracks for 1.5, 3, 6, 9 and 
15 M G (black full lines), 0.01, 0.1, 1 and 10 Myr isochrones (blue thin da- 
shed lines), and the ZAMS (black dot-dashe d line) are also plotted. The 
birthline taken from Behrend & Maeder 12001) is plotted with a blue thick 
dashed line. The convective/radiative phase transition is overplotted with an 
orange dot-dot-dot-dashed line. 

power of 65000. The data were reduced using the 'Libre Esprit' 
package especially developed for ES PaDOnS and Narv al, and ins- 
talled at the CFHT and at the TBL dDonati et alii 19971) . After re- 
duction we obtained the intensity (Stokes /) and the circular polari- 
zation (Stokes V) spectra of the stars observed. We then applied the 
least-squares deconvolution (LSD) procedure, a multi-line analysis 
technique, allowing us to increase the signal-to-noise ratio of our 
data by extracting the mean Stokes / and V profiles for each spec- 
trum. The Stokes / profiles were then used to measure the v sin i at 
the surface of each star. More details on the reduction and analysis 
procedures are given in paper I. 

Our sample contains 70 field Herbig Ae/BeQ. In paper I we 
used the photometric data and distances of the stars, compiled from 
the literature, to give an estimation of their luminosity. By fitting 
our observations with synthetic spectra, we measured their effec- 
tive temperature. We then placed the stars in an HR diagram, and 
compared their position with PMS evolutionary tracks computed 
using the CESAM code, version 2K dMorelll 19971) with solar abun- 
dances, in order to derive the mass, radius and age of each star, as 
well as the expected radius of the star when it reaches the ZAMS 
(Table 2 of paper I). For one star, HD 34282 , that has a low metalli- 
city compared to the sun dMerfn et alj2004h, we used the age, mass 
and radius determination of iMerin et alj d2004l) . however we were 
not able to estimate the ZAMS radius, because our CESAM models 
are of solar metallicity. Therefore we have rejected this star in the 
following analyses. More details on these procedures are given in 
Paper I. In Fig.Q]the stars are plotted in a HR diagram superimpo- 
sed on PMS evolutionary tracks and isochrones. The mass distri- 
bution of our sample, as well as the masses as a function of ages, 



few of them are now known to belong to associations such as Orion, but 
were thought to be field stars at the beginning of this study 
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Figure 2. left : Mass distribution of the sample, right : Masses of the stars 
as a function of their age. The red dashed line is a fit to the data points for 
masses lower than 5 M Q . 



are plotted in Fig. [2] For two stars (HD 98922 and IL Cep) that fall 
outside of the PMS region of the HR diagram, we couldn't estimate 
their mass, radius and age. Therefore they won't be included in the 
following analyses. 

In Fig. [T] we incl ude the birthline as computed by 
iBehre nd & Maeder ( 200ll. BM01 hereinafter). For reasons ex- 
plained in Paper I, we favoured their birthline to those of 
IPalla & StahleJ dl993h . The birthline of BM01 has been compu- 
ted by assuming a modulated accretion rate during the proto-stellar 
phase that increases as the star gains in mass. Once the close envi- 
ronment becomes poor in gas, the strong accretion phase stops and 
the star evolves on the PMS track from the birthline to the zero-age 
main-sequence (ZAMS). The age of each star has been computed 
from the birthline. This age indicates therefore the time that the 
star spent on the pre-main sequence phase, and does not include 
the time of the formation of the star. As a result some stars, situa- 
ted in the HR diagram very close to the birthline, can have ages as 
young as few thousands of years. 

We have plotted in Fig.[TJthe limit after which the mass of the 
convective envelope is less than 1% of the mass of the star. This 



limit represents the end of the convective phase and the start of 
the radiative phase of PMS evolution, during which the star is pre- 
dominantly radiative (a small convective core appears only at the 
very end of the PMS phase, when the luminosity starts to decrease 
again). We observe in Fig.Q]that the less massive stars of the sample 
are more concentrated toward the ZAMS and that the contrary i s 
observed at high mass, as already reported bv lVieira et al.l d2003l) . 
We suggest that this is due to a selection effect, using Herbig Ae/Be 
criteria. At high mass the close stellar environment dissipates fas- 
ter than at low mass, du e to strong winds and photoionisation from 
the stellar UV flux (e.g. iMeeus et al.ll200ll ; iDullemond et alfeOOll : 
lAlonso-Albi et alj2009h . Therefore, if stars are selected using Her- 
big Ae/Be criteria that probe a CS environment rich in dust and 
gas, we are more likely to select stars at earlier PMS evolutionary 
stages in massive stars than in lower mass stars, before the the CS 
environment has dissipated sufficiently to become undetectable in 
the IR or in the spectrum. We argue that PMS massive stars close 
to the ZAMS could have already dissipated their CS environment. 
As a result they do not show anymore Herbig Ae/Be criteria, and 
therefore cannot be part of our sample. 

In the case of low-mass stars, the whole sample is concentra- 
ted in the radiative phase of PMS evolution, as illustrated by the 
dearth of points on the right of the convective/radiative phase tran- 
sition in Fig. 1. This part of the HR diagram is in fact populated 
with the intermediate-mass T Tauri stars (IMTTS). The IMTTS are 
the evolutionary progenitors of th e Herbig Ae/Be s tars with spectral 
types ranging from K to F (e.g. lCalvet et al. i 2004), which explains 
why they have not been included in our sample. All these reasons 
explain why a strong correlation is found between the mass and the 
age of the stars of our sample. 

In Fig. 2 (bottom), the age-mass correlation is evident for 
masses lower than 5 M G . However, between 5 and 10 M G , while 
the stars are younger than at lower mass - which is expected as 
above 5 M G the PMS lifetime is lower than 0.4 Myr - the ages 
are highly scattered and no correlation of the ages at these masses 
is evident. In order to estimate the effect of the correlation found 
at lower mass on the analyses described in Section 3, we have 
performed a linear fit of the data points. We find a relation : 
M/M Q = 3.6 - 1.9 logOVlMyr), where t is the age. 



3 ANALYSIS OF THE PROJECTED ANGULAR 
MOMENTA OF THE SAMPLE 

One of the aims of this survey is to understand angular mo- 
mentum evolution during the PMS phase of intermediate mass 
stars, and to estimate the role of the magnetic field in this evolu- 
tion. We present below two different analyses to address these ques- 
tions. The first aims to find a correlation of the angular momentum 
of a sample of HAeBe stars with their age, while the second com- 
pares the v sin i distributions of the same HAeBe sample to various 
samples of MS A/B stars in order to investigate if angular momen- 
tum evolution is expected in the future of HAeBe stars, before they 
reach the ZAMS (Section 4). 

Today, we have no knowledge of the internal rotation profile 
of Herbig Ae/Be stars. During its evolution, the rotation profile of 
a PMS star is expected to change as the star contracts and as angu- 
lar momentum is transported through various mechanisms not fully 
understood yet. While understanding the behaviour of internal ro- 
tation in PMS stars is an important field of theoretical research, it is 
also becoming observationally accessible thanks to the CoRoT and 
Kepler satellites. However, our knowledge of internal rotation on 
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the PMS is still very poor. We therefore can consider two extreme 
hypotheses : (i) a solid body rotation (SR), and (ii) a constant speci- 
fic (i.e. per unit of mass) angular momentum (CS), between which 
the real internal rotation profile would be. The solid body rotation 
(hypothesis i) assumes that the transport of angular momentum in- 
side the star is highly efficient all along the PMS evolution. This 
hypothesis can be considered as reasonable inside a fully convec- 
tive star, such as low-mass stars at the beginning of the PMS phase. 
However at intermediate-mass, during the major part of their PMS 
evolution, the stars are mainly radiative. Hypothesis (?) is therefore 
too far from the reality for Herbig Ae/Be stars. On the contrary, if 
we consider an isolated star without exchange (of matter or angu- 
lar momentum) with its environment and without internal transport 
of angular momentum, we can assume a constant specific angular 
momentum. As the star contracts, the mass will get more and more 
concentrated into the core and the angular momentum will tend to 
behave as 1 /r 2 . In a radiative interior we would expect the transport 
of angular momentum to be much less efficient than in a convective 
interior and therefore hypothesis (ii) would be more reasonable for 
the radiative phase of the PMS evolution. However, in reality du- 
ring the PMS phase a star will exchange matter, and therefore angu- 
lar momentum with its environment. Internal transport of angular 
momentum through various processes, such as pulsations or meri- 
dional circulation, are also expected to be active with an efficiency 
that is unknown. The real internal profile is therefore expected to 
lie between both hypotheses. 

We have calculated the angular momentum of each star in our 
sample using the internal density profile at the previously estima- 
ted age, using CES AM models, and both assumptions on the stellar 
rotation profile. In the case of a spherically symmetric star, the an- 
gular momentum element dj is given by : 
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dj = —r'iL)(r)dm, 



(1) 



where dm is the mass of a spherical shell inside the star at a distance 
r from the center, and u>(r) is the angular velocity of the shell. 

In the case of solid body rotation, the angular velocity is 
constant as a function of r and equal to the angular velocity at the 
surface of the star : co(f) = Q.,, Vr < R,. The total stellar angu- 
lar momentum is therefore computed by integrating the following 
equation throughout the volume of the star : 



Jsr - 



yQ» j p(r)r 4 d r , 



(2) 



where p(r) is the density of the star as a function of radius. 

In the second case, the specific angular momentum dj/dm is 
constant and equal to \R 2 ,Cl, at the stellar surface. The angular ve- 
locity is therefore dependent on r and equal to : 



IR \ 2 
co(r) = J-i) a, 

and the angular momentum is computed as follows : 
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RlCl, p{r)?dr. 



(3) 



(4) 



The angular velocity at the surface of the star is computed 
using the estimated radius R, and the projected rotational velocity 
v sin i. As sin i is unknown, the angular velocities, hence the angular 
momenta, of the stars are computed up to a factor of sin i. We as- 
sume that whatever the ages, masses, or situation within the Galaxy 
of large samples of stars, the distributions of sin i can be considered 
similar in all cases. We therefore expect that if all quantities consi- 



Figure 3. Angular momentum of the sample as a function of age, assuming 
a solid body rotation with a normal y-axis scale (upper) and a logarith- 
mic y-axis scale (bottom). The circles represent the normal stars, while the 
squares represent the magnetic stars. The symbols are defined as follows : 
large symbol : M > 10 M , medium filled symbol : 7 < M < 10 M Q , 
medium open symbol : 5 < M < 1 M Q , small symbol : M < 5 M Q . The 
black lines plotted in the bottom panel are the predicted /sr of a sample of 
stars evolving at constant angular momentum, with a mass-age correlation 
similar to that of our sample (see Section 2), and with vsini of 10 km.s~' 
(dashed), 50 km.s~' (dot-dot-dot-dashed), 100 km.s~' (long-dashed), and 
250 km.s~'(dot-dashed). The bumps observed around 1 Myr mark the tran- 
sition between fully radiative stars (on the left of the bump), and stars with a 
convective core (on the right). The upper or lower error bars in J that appear 
very small, are for stars very close to the birthline or the ZAMS. 



dered in one statistical analysis are computed or measured up to a 
factor of sini, it won't affect the result. 

We have divided our sample into two distinct sub-samples, 
one containing the magnetic HAeBe stars and the other containing 
the non-magnetic, non-binary stars (or "normal stars"). The magne- 
tic sample contains the five confirmed magnetic HAeBe stars (HD 
200775, V380 Ori, HD 190073, HD 72106 and LP Ori) that have 
been reported in paper I. The normal sample contains the 58 stars 
in which no magnetic fields have been detected and which are not 
member of a close binary system (paper I). From this sample we 
exclude MWC 1080 that has no vsini HD 250550 that has no es- 
timation of its fundamental parameters, and HD 34282, that has no 
estimation of a ZAMS radius. The final normal sample considered 
in this study contains 55 stars. Both magnetic and normal samples 
will be compared, then treated separately, as the magnetic fields 
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are expected to have a strong influence on the rotational evolution 
of the stars. The non-magnetic binary stars have not been included 
into the analysis as tidal interaction between both components is 
expected to modify the stellar angular momentum. 



3.1 The magnetic stars 

The favoured hypothesis explaining the slow rotation of ma- 
gnetic Ap/Bp stars is magnetic braking occurring during the PMS. 
This would slow magnetic stars much more than non-magnetic 
stars, throug h the interactio n of the stellar magnetic field with its 
environment dStepiehl200rl) . According to this hypothesis a change 
of the angular momentum of PMS A/B stars is expected between 
magnetic and non-magnetic stars. To test this hypothesis we have 
plotted, in Fig. [3] the solid body rotation angular momentum (Jsr) 
as a function of age for the magnetic (squares) and normal (open 
and filled circles) HAeBe stars. The size of the symbol is indicative 
of the stellar mass. We don't observe any correlation of the angular 
momentum of the magnetic stars with age. The same is observed 
when a constant specific angular momentum is assumed instead. 
However, when considering the small number of magnetic stars in 
this sample, the range of stellar mass covered by this sample, and 
the fact that the angular momenta have been calculated up to a fac- 
tor sin i, if a correlation exists, it would be difficult to detect. 

On the other hand, even if the sample is small, it is worth no- 
ting that four out of the five magnetic stars have an angular momen- 
tum much lower than the normal sample. The same is observed in 
the case of a constant specific angular momentum. 



3.2 The non-magnetic stars 

Even in the absence of strong magnetic fields, HAeBe stars 
can experience an exchange of matter with their environment 
through accretion disks and/or winds, which are observable through 
the emission properties of their spectra. Therefore the stars can ex- 
change angular momentum with their environment. Furthermore, 
while the magnetic stars are those with strong detected magnetic 
fields, there is still the possibility that, in the sample of the "non- 
magnetic" stars, some stars could host magnetic fields too faint to 
be detected (see Paper III). If such fields exist they could enhance 
the loss of angular momentum through an exchange of matter oc- 
curring during the PMS phase of their evolution. We therefore per- 
formed an analysis of the angular momentum of the non-magnetic 
stars in order to understand their angular momentum evolution du- 
ring the PMS phase. 

In Fig.[3] while the angular momenta of the normal stars gene- 
rally spread from J$ R = to a maximum value, which is due to the 
uncorrected sin i factor, the maximum value clearly decreases with 
age whatever the mass. In the non logarithmic plot (Fig. 5 top), this 
trend is most significant in the middle sub-sample with masses bet- 
ween 5 and 10 M Q (medium open and filled circles). We checked 
that at these masses the trend is not due to the age-mass correlation 
pointed out in Section 2 by observing that the graph area covered 
by a 5-7 M G sub-sample (open medium circles) is approximately 
the same as the area covered by a 7-10 M Q sub-sample (filled me- 
dium circle). Furthermore, between 5 and 10 M G , as pointed out in 
Section 2, no correlation of the age and the mass is observed. 

In the logarithmic plot of the angular momentum (Fig. 3 bot- 
tom), a clear decrease of the angular momentum of the stars with 
masses below 5 M Q (small open circles) is observed. However 
at these masses, a strong age-mass correlation is present in our 
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Figure 4. vsini histograms of the magnetic (left) and the non-magnetic 
HAeBe stars. 



sample. In order to check if it can explain this decrease, we have 
calculated the expected angular momentum of a simulated sample 
of stars of various masses between 1.5 and 5 M Q , evolving at 
constant angular momentum, with ages following the correlation 
found in Sec. 2, and for four vsini : 10, 50, 100 and 250 km.s -1 . 
We observe (Fig. 3 bottom) that the trend is well reproduced by the 
angular momentum of the simulated sample, and that most of the 
lower mass stars lies between the curves of 50 and 250 km.s~' in 
v sin i, values that are consistent with the v sin i of most of our stars. 
Therefore, the age-mass correlation of our sample at lower mass 
can explain alone the decrease of the angular momentum observed 
in Fig. 3 (bottom). 

In order to test the impact of the choice of the internal rota- 
tion profile, we did the same plots using the total stellar angular 
momentum computed by assuming a constant specific angular mo- 
mentum. While the total angular momentum of the stars is gene- 
rally found to be larger than in the case of a solid body rotation, 
we obtain a similar plot, with the same general trend of the maxi- 
mum values decreasing with time for stellar masses above 5 M Q . 
At lower mass, a trend similar to the SR case is observed and is 
also fully explained with the age-mass correlation that our sample 
shows at these masses. We therefore find that, whatever the stellar 
internal rotation profile considered, the highest mass (M > 5 M Q ) 
PMS stars seem to loose angular momentum as they evolve towards 
the ZAMS, while the lower mass stars seem to evolve with constant 
angular momentum. 



4 ANALYSIS OF THE PROJECTED ROTATIONAL 
VELOCITIES 

4.1 The magnetic stars 

A strong dichotomy is observed in t he rotation velocities of 
non-magnetic and magnetic MS A/B (e.g. lAbt&Morreiill 19951) : 
the former are fast rotators, while the latter are much slower ro- 
tators (with rotational periods generally longer than one day). In 
order to understand this dichotomy and its evolution, we have plot- 
ted, in Fig. [4] the v sin i distributions of the magnetic (right) and 
non-magnetic (left) HAeBe stars. We observe that already during 
the PMS phase there is a dichotomy similar to that observed among 
MS A/B stars. All magnetic HAeBe stars rotate with v sin i lower 
than 100 km.s -1 . Since the number of magnetic stars is small, one 
could argue that the low v sin i seen on these stars could be due to a 
low inclination of the rotation axis. However, the magnetic analysis 
of the individual stars allows for a direct determination of their ro- 
tation periods, independent of the inclination. We find periods ran- 
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ging from 1.6 to 4.3 d ays jAlecian et alj2008 d ; lFolsom et alj2008l : 
Alec ian et ai] |2009b). We can therefore conclude confidently that 
all magnetic HAeBe stars of our sample are slow rotators, while 
only part of the non-magnetic stars of our sample rotate slowly (or 
have a low v sin i). 

We performed a p- value statistical test to check the hypothesis 
that the v sin i of the magnetic sample can exist by chance, if we 
assume that the v sin i of the magnetic and normal sample follow 
the same distribution. In this aim we have comp uted the Student' s 
f-value and the associated significance p-value dPress et al.|[l992h . 
We find a p-value of 0.002, that is much lower than 0.05, mea- 
ning that the null hypothesis - that the observed v sin i values of the 
magnetic stars can be ascribed to chance alone - is rejected. The 
differences between both distributions are therefore real and have a 
physical meaning. 

4.2 The non-magnetic stars 

The rotational v elocities of MS A/B stars hav e already been 
largely measured (see Rov er. Zorec & Gomez 2007, and references 
therein), and their global distribution is known. We can therefore 
search for a potential evolution in the angular momentum during 
the PMS by predicting the rotation velocities that HAeBe stars will 
have once they will reach the ZAMS and comparing them to the 
distribution of MS v sin i. There are many reasons to believe that 
the evolution of angular mome ntum in A/B stars is negligible du- 
ring the MS (e.g. Stepieri 2000). Furthermore, as the evolution of a 
MS star is slower during the first half of the MS, most MS stars are 
closer to the ZAMS than the terminal age main sequence (TAMS). 
Hence, it is reasonable to assume that the stars from the MS ca- 
talogues discussed below can be treated approximately as ZAMS 
stars. The ZAMS v sin i, that we have computed, can therefore be 
directly compared to the v sin i of general MS stars. 

Assuming a uniform angular momentum evolution through 
the PMS, we have estimated the rotational velocities on the ZAMS 
for our sample as follows : 

. . . .^zams I* 
vzams sin ; = v sin i — — (5) 

'ZAMS 

where, I, and /zams are the moment of inertia at the current stellar 
age, and on the ZAMS, respectively. The moments of inertia have 
been computed by integrating the following equations for both stel- 
lar rotation profile assumptions : 

/sr = y J P(r)r 4 dr (6) 

Ics = Y Rl * J P^Vdr (7) 

For both angular momentum treatments, the resulting distributions 
of ZAMS i' sin ( show clear differences when compared to the PMS 
vsku one : they are shifted towards larger vsin; and are more 
spread. It is the result of a decreasing of the radius and an increase 
of the density. 

We have considered different samples of rotational velocities 
of MS stars to compare our results with. The following describes in 
detail the comparisons with both angular momentum treatments. 

4.2.1 Comparison with the CPRV catalogue 

We us ed the Catalogue of Pro j ected Rotational Velocities 
(CPRV) of iGlebocki & Stawikowskl (2000) which is a compila- 
tion of the observed v sin i of over 1 1000 stars of all spectral types 



and luminosity classes. This catalogue includes numerous classes 
of stars that should not be taken into account in our analysis for 
reasons explained below. We have selected the suitable stars as fol- 
lows : 

(i) we excluded all stars with uncertain v sin i (as indicated in 
the catalogue) : the number of stars in the catalogue is sufficiently 
large to keep only well constrained v sin i. 

(ii) when multiple v sin ( have been published for a star, we kept 
the more accurate value (i.e. the one with the smallest error bare). 

(iii) we only kept stars with luminosity class V which can be 
compared to the ZAMS v sin i of HAeBe stars. 

(iv) we excluded binary systems, as identified in the CPRV : 
tidal interactions between the components of binary systems lead 
to a strong rotational braking. These stars should therefore not be 
taken into account in our test of angular momentum conservation. 

(v) we exclude all peculiar stars : all Ap/Bp stars are magnetic 
and should be rejected, as we are considering only non-magnetic 
stars in this section. Among the other class of peculiar stars, Am 
stars exhibit a binary incidence close to 100% ; for reasons explai- 
ned above, we exclude them too. 

(vi) we only kept intermediate mass stars : in order to deter- 
mine the masses of MS st ars, we used the Stellar Mass Catalogue 
(SMaC) o f lBelikovl l fl995T) . that contains the dynamical determina- 
tion of the masses of binary systems of luminosity class V, as well 
as their spectral types. By plotting log(M/M Q ) as a function of spec- 
tral type, we found two linear regressions, one for spectral types 
earlier than AO and one for spectral types later than AO (cf. Ap- 
pendix A). Using these calibrations and the tabulated spectral type 
in the CPRV, we determined the masses of the stars. Only masses 
between 1.5 and 20 M Q have been used. 

Some binary systems are indicated in the CPRV, but in order 
to exclude all of them, we used the most rece nt version of t he Wa - 
shington Double Stars (WDS) catalogue of iMason et afl l l200ll) . 
We first rejected from the WDS catalogue, all dubious binaries or 
multiple systems and all non-physical systems, as indicated in the 
catalogue itself. Then using the coordinates of the CPRV and the 
WDS catalogues, we calculated the distances between the stars of 
both catalogues. Whenever the distance was lower than 90 arcsec 
(which is roughly the precision on the coordinates of the CPRV 
stars) the star was considered as double and was rejected. The final 
number of stars of the CPRV suitable for our analysis is 1823 stars. 

The mass range of our HAeBe sample is broad (from 1.5 to 
17 M Q , see Table 2 of Paper I). Many phenomena during the PMS 
phase are strongly sensitive to the mass (e.g. stellar winds, accre- 
tion disks, PMS lifetime). Therefore we divided the ZAMS HAeBe 
and MS samples into three sub-samples of different mass ranges, in 
order to keep the number of stars equal in each bin : 

- #1 : M < 2.3 M , 

- #2 : 2.3 < M < 3.5 M , and 

- #3 : M > 3.5 M Q . 

Each one of the ZAMS HAeBe sub-samples contains 19, 19 and 
17 stars respectively, while the MS CPRV sub-samples contain 
853, 415 and 567 stars, respectively. We first compared the 
v sin i distributions of the ZAMS HAeBe and the MS stars for the 
three mass ranges. However, while the two first sub-samples cover 
a relatively low-mass range (compared to the typical error bars, see 

2 The catalogue has been downloaded from the WDS website : 
http ://ad. usno.navy.mil/wds/ and the last update was Oct. 2 nd 2008 
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Figure 5. W sin i histograms of the HAeBe stars projected on the ZAMS (blue field bares) assuming a solid body rotation (top) or as suming a constant specific 
angular momentum (bottom), and of the MS stars from the catalogue CPRV (red empty bares) (Glebocki & Stawikowski 2000), for three different mass 
ranges : M < 2.3 M Q , 2.3 < M < 3.5 M Q , M > 3.5 M G . 



Fig. 1), the third sub-sample covers a wide mass range that need 
to be taken into account when comparing two v sin i distributions. 
The mass distributions of the HAeBe and CPRV samples are dif- 
ferent : a lack of high-mass stars (above 9 M Q ) and an excess of lo- 
wer mass stars (below 5 M Q ) is observed in the HAeBe distribution 
with respect to the CPRV one. Due to these differences, comparing 
two v sin ! distributions is meaningless as the potential differences 
could be due to either a difference in mass (or radius) or in an- 
gular momentum. Therefore, in order to minimise this effect, we 
have chosen to compare the rotation rates (Wsini), instead of the 
v sin i, which are defined as the ratios of the v sin (' to the break-up 
velocity : 

v sin i 

W sin / = — (8) 

(GM/ i?ZAMS ) 

where M and R are the mass and radius of the star, and G is the gra- 
vitational constant. The masses of the HAeBe sample have been de- 
termined in paper I. The mass of the CPRV sample have been esti- 
mated by using a mass-radius relation. This relation has been deter- 
mined by using the Catalogue o f orbital elements, mass and lumi - 
nosities of close double stars of ISvechnikov & Bessonovi] j 1984b . 
In this aim we first performed a selection on the stellar parameters 
as follows : 

- we kept only stars of luminosity class V 

- we excluded all stars with uncertain or unknown mass ratio, 
mass, or radius 

- we kept stars more masive than 1 M . 

Then we plotted the log(R/R B ) as a function of log(M/M Q ) and 
performed a linear regression (see Appendix B). 

In Fig. [5] are plotted the distributions of the rotation rates 
W sin i of the ZAMS HAeBe and the MS CPRV stars for all masses, 
for the three mass ranges, and for both angular momentum as- 
sumptions (SR and CS). Whatever the mass, the ZAMS SR and 
the CPRV look similar. However, when comparing the ZAMS CS 
and CPRV distributions, all masses together, we note some diffe- 
rences : the first one looks slightly shifted towards larger W sin i, 
and it shows a lack of slow rotators in favour of fast rotators. These 
differences seem to be reproduced in the third mass sub-sample. 

In order to determine if the differences that our eye seem to 
note are significant, we performed two statistical tests. Both tests 



estimate the significance level of the rejection of the null hypo- 
thesis (i.e. both distributions are the same). Fir st, we performed a 
Kolmogorov-Smirnov (KS) test as described bv lPress et aT1jl992h . 
This test compares the cumulative distributions of both samples. It 
is sensitive to the shift of one distribution with respect to the other. 
When we apply it to each panel of Fig. 5, the significance levels 
are all below 95 % except for the panels g (99%) and h (98%). It 
confirms that real differences exist between the ZAMS CS and the 
CPRV catalogue at high mass (M > 3.5M Q ). In the case of a solid 
body rotation, the highest significance level is at high-mass, with 
a value of 93%. However we consider that 93% is not statistically 
significant (typical values would be between 95 and 99%). 

The second test t hat we performed is a chi-square test, as des- 
cribed bv lPress et ail dl992h . This test compares two distributions 
bin per bin by calculating the chi-square between both distributions. 
By performing a standard chi-square test we estimate the probabi- 
lity of rejection of the null hypothesis. As with the KS test, we find 
that the ZAMS HAeBe distributions in the panel h only is signi- 
ficantly different to the CPRV distribution at a level of 99.993%. 
In the case of a solid rotation, the highest probability of rejection 
is found at high mass (panel c) with a value of 87%, which is not 
enough to confidently reject the null hypothesis. 

Both tests find that the differences between the ZAMS CS 
and the CPRV catalogue are significant at masses above 3.5 M Q , 
indicating that we predict too much fast rotators by assuming a 
constant angular momentum. In the case of solid body rotation, 
above 3.5 M Q , it also seems that an excess of fast rotators is predic- 
ted, compared to MS stars. However neither statistical test is able 
to quantitatively confirm this 



4.2.2 Comparison with the RZG catalogue 

In order to verify that our above conclusions are not sample- 
dependent, we performed the same test with the catalogue of rota- 
tional velocities of A-type stars of iRover. Zorec & Gomej J2007L 
RZG hereinafter). We rejected all the binaries and chemically pe- 
culiar stars from this sample, as indicated in the catalogue itself, but 
also as indicated in the WDS catalogue. The final number of stars 
from this catalogue, considered in this study, is 981. The mass of 
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Figure 6. W sin i histograms of the HAeBe projected on the ZAMS (blue 
field bares) assuming a solid body rotation (top) or assuming a constant spe- 
cific angular momentum (bott om), and of the MS stars from the catalogue 
RZG (red empty bares) iRover et al.ll2007l) . for two different mass ranges : 
M < 2.15 M , 2.15 < M < 2.6 M Q . 




c) 



0.00.20.40.60.81.0 
Wsini 



d) 



0.00.20.40.60.81.0 
Wsini 



Figure 7. Wsini histograms of the HAeBe projected on the ZAMS (blue 
field bares) assuming a solid body rotation (top) or assuming a constant spe- 
cific angular momentum ( bottom), and of the MS stars from the catalogue 
ffPer (red empty bares) of lProsse 3 il992h . for two different mass ranges : 
M < 2.3 Mq, 2.3 < M < 3.5 M G . 



the stars in the RZG catalogue are all below 2.6 M Q . We therefore 
divided both ZAMS HAeBe and RZG samples in two sub-samples : 
M < 2.15 M Q and 2.15 < M < 2.6 M Q . The numbers of stars are 
13 in each of the ZAMS HAeBe sub-samples, and they are 496 
and 485 in the RZG sub-samples, respectively. The distributions 
are plotted in Fig. [6] In spite of some irregular ZAMS HAeBe dis- 
tributions, we don't detect with the eye any systematic differences 
between the distributions. The KS test gives significance levels be- 
low 95% for all panels. Therefore, whatever the treatment of the in- 
ternal angular momentum, no difference is found between the RZG 
and ZAMS HAeBe samples. 



perform a reliable statistical analysis. We therefore only consider 
the two first mass ranges : M < 2.3 M G and 2.3 < M < 3.5 M G . 
The distributions are plotted in Fig. [7] The ZAMS and MS distri- 
butions are highly irregular, imputable to small samples. In spite of 
that panel b) seems to show some differences : a lack of fast rota- 
tors and an excess of slow rotators is predicted in the ZAMS SR 
sample. The KS test gives significance levels below 95 % in all pa- 
nels except in panel b, where it finds 96 %, in contradiction with 
the results obtained with the CPRV and RZG samples. 



4.2.3 Comparison with the a Per open cluster 

Finally, we performed a third comparison with the W sin i dis- 
tributions of the inter mediate mass st ars in the open cluster a Per, 
from the catalogue of IProssetl |l992). These stars are co-evolved 
and the age of the cluster is estimated to be around 80 Myr. The 
duration of the PMS phase for intermediate mass stars is 30 Myr 
maximum (see Table 2 of Paper I). The a Per sample is there- 
fore more representative of intermediate mass ZAMS stars than the 
other samples used above. From the original a Per catalogue we 
performed the following selections : 

(i) we excluded all stars with uncertain v sin i. 

(ii) we excluded binary systems, using the WDS catalogue. 

(iii) we only kept intermediate mass stars : we used the same 
method as the CPRV catalogue to determine the masses. This me- 
thod requires a kn own spectral typ e. If the spectral type is unknown 
in the catalogue of lProssen il992h , we used the one indicated in the 
SIMBAD Astronomical database[| 

We divided both the ZAMS HAeBe and a Per samples into 3 
sub-samples of the same mass ranges as the CPRV ones. The num- 
bers of stars in the a Per sub-samples are 30, 20, and 9, respectively. 
Below 10 stars, we consider that the number of stars is too low to 

3 http ://simbad. u-strasbg.fr/simbad/ 



5 CONCLUSIONS AND DISCUSSION 

This paper is the second in a series that presents the results 
of a high-resolution spectropolarimetric analysis of a sample of 70 
Herbig Ae/Be stars. While paper I describes the data and paper III 
will examine the magnetism of Herbig Ae/Be stars, here we address 
the problem of stellar angular momentum evolution during the pre- 
main sequence phase at intermediate mass. 

Using the vsin; values published in paper I, we performed 
two separate analyses. In both analyses we have used two internal 
rotation profile assumptions : a solid rotation, and a constant speci- 
fic angular momentum. As the internal rotation profiles of Herbig 
Ae/Be stars is unknown, we chose these two extreme hypotheses 
for our analysis. In the first analysis we have computed the angular 
momentum of each star of our sample to search for a correlation 
with the age for the magnetic sample on one side, and the normal 
(i.e. non-magnetic and non-binary) sample on the other side. In the 
second analysis, we first compared the vsini distributions of the 
magnetic to the normal samples, then we predicted the v sin i that 
our normal sample will have on the ZAMS using both above as- 
sumptions, and compared the derived distributions to published 
v sin ( distributions of intermediate mass MS stars. In both analyses 
we have analysed the normal sample within mass ranges, in order to 
avoid mass-dependent conclusions. However, in order to minimise 
the differences found between both ZAMS and MS mass distribu- 
tions above 3.4 M Q , we have derived the rotation rates W sin i (i.e. 
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the v sin i over the escape velocity) of the samples, and compared 
the W sin i distributions (instead of the v sin (). 

5.1 Magnetic field impact on the angular momentum 
evolution 

Both analyses revealed that the magnetic Herbig Ae/Be stars 
of our sample are much slower rotators than the normal ones, as 
observed on the main sequence, suggesting that the magnetic stars 
have experienced more braking than the normal stars. Table 2 of 
paper I provides the ages of the stars and the duration of their PMS 
phases. We observe that, while some of the magnetic stars have 
completed a large fraction of their PMS phase, one of them (HD 
190073) has only completed a third of its PMS phase. Furthermore, 
HD 200775 has only an age of 16000 years. These stars are very 
young and are already rotating very slowly. If magnetic braking 
is responsible for these slow rotators, it must act very early, and 
relatively fast (in less than 16000 years), during the PMS phase, or 
even before the PMS phase. 

We could propose an alternative to the magnetic braking 
theory. We observe, in a large number of H AeBe stars, wind signa- 
tures in their UV and optical spectra (e.g. Finkenzeller & Mundt 
1984, paper III). Winds seem therefore frequent during the PMS 
phase at intermediate mass. A very young star, newly formed, and 
starting its PMS phase is rotating very fast and experiences mass 
loss, and therefore angular momentum loss, due to its wind. The 
surface of the star is therefore rotating more slowly than the core. 
A velocity gradient appe ars between the surfac e and the core, indu- 
cing sheared turbulence dLignieres etal]|l996h . We argue that this 
turbulence could be sufficiently strong - depending on the rotational 
energy compared to the magn etic energy - to dissip ate the magnetic 
fields inside very young stars dAuriere et al .120071) . Therefore, only 
relatively slow rotators could keep their magnetic fields. However 
this theory requires that a population of slowly rotating stars exists 
before they start their PMS evolution. 

Our analysis does not allow us to favour either of these hy- 
potheses for the following reasons. First the number of known ma- 
gnetic HAeBe stars is very small. It is therefore difficult to per- 
form any statistical analysis that could reveal a dependance of their 
angular momentum with magnetic strength, age, environment, ini- 
tial conditions, or stellar properties (e.g. masses or rotation per- 
iod). Secondly, we need to better characterise the time-scales of the 
magnetic braking, or magnetic dissipation which could help us to 
favour one of the hypotheses, when compared to the ages of the 
stars. Third, the uncertainties on ages, initial conditions and envi- 
ronment are very large among field HAeBe stars. It is necessary to 
perform the same analysis on HAeBe star members of very young 
clusters of different well determined ages, in which sets of stars 
have the same age and same initial conditions. Such a study could 
help us to disentangle evolutionary effects from effects of the ini- 
tial conditions. We have collected a series of spectropolarimetric 
data of ~ 60 HAeBe stars in three very young clusters (NGC 2244, 
NGC 2264 and NGC 6611), in w hich we discovered 3 magnetic 
stars dAlecian et alj|2008bl.l2009al. Alecian et al. in prep.). An ana- 
lysis of the whole sample is in progress and will be published in the 
near future. 

Finally three of the five magnetic stars constituing this sample 
belong to binary systems, therefore tidal braking could also be res- 
ponsible of their slow rotation. However the three systems are not 
synchronised, and are wide (a/R > 10, where a is the semi- major 
axis o f the orbit, and r is the radius of the star). According to lZahrj 
(1977), in the case of wide binaries, the evolution of the angular 



velocity tends towards an evolution in absence of tidal braking, 
meaning that tidal braking in wide binaries are negligible, espe- 
cially during the first stages of stellar evolution. It implies that 
in such wide binaries, the synchronisation would occur on times- 
cales of the order of the main-sequence lifetime, consistent with 
the observation of our three very young and not-synchronised sys- 
tems. We therefore think that it is unlikely that tidal torques have 
been efficient enough to brake these young stars. Furthermore, in 
these systems, the magnetic component is systematically rotating 
slower than its non-magnetic companion. One of both hypotheses, 
involving magnetic field, that we exposed above, appear therefore 
more reasonable to explain the slow rotation of the magnetic Her- 
big Ae/Be stars. 



5.2 Angular momentum evolution of the normal sample 

First, we find that whatever the choice of the internal angular 
momentum treatment (solid body rotation or constant specific an- 
gular momentum), the results are the same. The plots of the angular 
momentum as a function of age show a clear trend at all masses. We 
checked that below 5 M G this trend can be fully explained with the 
age-mass correlation that our sample shows at lower mass. On the 
contrary, we find that above 5 M Q , the observed trend is real : the 
angular momentum of the stars decreases as the age increases. 

We have compared the rotation rates Wsin; of our ZAMS 
samples with thos e of the normal MS intermediate -mass stars in the 
CPRV catalogue dGlebocki & Stawikowskl 2000). We find that in 
the case of a solid body rotation, whatever the mass, no significant 
differences are observed between ZAMS HAeBe and MS CPRV 
samples. This result does not allow us to confirm the decrease of 
the angular momentum with age described above, but neither is it in 
contradiction with this observation (a KS- test can only reject the 
null hypothesis, it cannot confirm it).. Furthermore we note with 
the eye an excess of fast rotators in our ZAMS sample at high-mass 
compared to the MS sample. 

In the case of a constant specific angular momentum, when all 
masses are considered, it looks like too many fast HAeBe rotators 
are predicted. This trend seems to be observed only at high masses 
(M > 3.4 M Q ), which is confirmed with a KS test. The comparison 
of our ZAMS HAeBe sample with two other MS samples (RZG 
and a Per) confirms the trend at low-mass (M < 3.4 M Q ). Unfor- 
tunately, the RZG and oTer have no or an insufficient number of 
high-mass stars (M > 3.4M Q ) to test if our conclusions are depen- 
dant on the choice of the MS sample. 

In conclusion, both analyses suggest that, whatever the in- 
ternal angular momentum treatment, the lowest mass stars (M < 
5 Mo) evolve with constant angular momentum during the PMS 
phase, while at higher mass, the stars lose angular momentum. 
A larger sample of stars, especially at high-mass, is required to 
confirm this result. 

A similar statisti cal analysis was performed in the past by 
iBohm & CataQ dl995l BC95 hereinafter) on a sample of 29 Herbig 
Ae/Be stars. Using MUSICOS (the ancestor of Narval and ESPa- 
DOnS) echelle spectra with a resolution of R = 38000, they mea- 
sured the v sin i of all the stars using two methods. The first method 
is fitting many spectral lines with a convolution of a rotation profile 
and Gaussian of instrumental width. The second method is a com- 
parison between many observed spectral lines and synthetic spectra 
of different rotational broadenings. They compared their v sin i to 
those of A and B stars in young open clusters. 
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Figure 8. v sin i of this work compared to the v sin i of BC95, for three mass 
ranges : M < 2.6 M Q (squares), 2.6 < M < 4 M Q (triangles), M > 4 M Q 
(dots). The dashed line is the first bisector. 



They concluded that, if the star rotates as a solid body, low- 
mass HAeBe stars (M < 2.6 M Q ) lose angular momentum, inter- 
mediate mass HAeBe stars (2.6 < M < 4 M Q ) evolve towards 
the ZAMS at constant angular momentum, while the high-mass 
HAeBe stars (M > 4 M G ) gain angular momentum. Except for 
the intermediate-mass stars, these results are in disagreement with 
oursQ 

By comparing our sample with that of BC95 we find 15 
stars in common. In order to understand the discrepancy bet- 
ween our conclusions and the conclusions of BC95, we compa- 
red our vsini measurements to those of BC95. In Figure [8] are 
plotted our vsin; as a function of BC95's. We observe that the 
v sin i measurements of the high mass stars from BC95 tend to be 
lower than ours, which could be the reason for the discrepancy. The 
reason of the differences obtained between BC95's vsin; and ours 
is not clear. It could come from the different treatment of the CS 
contamination. However independent v sin i determination from the 
same set of data t han ours gives v sin i values consistent with ours 
jFolsom et al]|2012l). 

IWolff et alj (2004), from their studies of PMS stars in Orion, 
concluded that intermediate-mass (1.5 < M < 3.5 M Q ) stars lose 
angular momentum before they start the PMS radiative phase, then, 
during the PMS radiative evolution they evolve at constant angular 
momentum. In Fig.[T]we observe that all the stars of our sample are 
experiencing the r adiative phase of th eir PMS evolution. Therefore, 
our work confirms Wolff et alj d2004f) 's conclusions. 

Our work suggests that the more massive normal HAeBe stars 
(M < 5 Mo) are losing angular momentum as they evolve towards 
the ZAMS. These stars have no strong magnetic fields, and they 
are most likely not members of close binary systems. In order to 
understand this loss, and also the observed difference between high 
and low-mass HAeBe stars, we have searched for signs of winds 



4 We checked that, if we take the same mass limits as BC95, our conclu- 
sions remain unchanged. 



in our spectra, especially in the Balmer lines, in the He I 5875 
A line, and in the multiplet 42 of Fe II. We also searched for si- 
gns of winds in UV spectra using the MAST archive, either in the 
C IV 1548/1550 A lines in massive stars or in the Mg II h & k 
doublet in the HAe stars, that can display strong signs of winds. Of 
the 11 stars above 5 M Q , we find that 8 of them show either clear 
P Cygni profiles or a strong blueshifted absorption component in 
one or many of the searched lines, revealing the presence of a wind 
in the close stellar environment. Winds seem to be largely present 
among massive HAeBe stars and could be at the origin of their loss 
of angular momentum. 

We have argued before that winds could also be partly respon- 
sible for the loss of angular momentum in magnetic stars. However, 
when coupled to strong magnetic fields, winds are expected to be 
much more efficient at carrying off angular momentum than in stars 
hosting faint or no magnetic fields. Therefore, even if all HAeBe 
stars have winds and are losing angular momentum, we still ex- 
pect magnetic stars to rotate more slowly than non-magnetic stars, 
which is consistent with our observations. . 

In the case of low-mass HAeBe stars, the predicted ZAMS 
vsini values are similar to those observed in MS stars. However 
we should note that most of the sample stars at low mass have 
completed more than 50% of the PMS phase. We therefore can- 
not conclude anything about the evolution of angular momentum 
during the first half of the PMS, for low-mass HAeBe stars. There 
are different ways to fill this gap : one way is to study HAeBe stars 
in very young clusters, such as NGC 6611, NGC 2244 and NGC 
2264, that have ages from less than 1 Myr to ~ 2.5 Myr. The results 
of our analysis of the spectropolarimetric data that we obtained for 
these clusters will be published in a future paper, and will discuss 
that point. Another way to fill this gap is to study intermediate-mass 
T Tauri stars, which have masses between 1.5 and ~ 3.5 M s , and 
are evolving through the convective phase of their PMS evolution. 
These stars have not been included into our analysis. An ESPa- 
DOnS program is in progress in order to investigate these stars, and 
complete our knowledge of PMS angular momentum evolution at 
intermediate mass. 
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Figure Al. Masses as a f unction of spectral types for the stars of the SMaC 
catalogue (Belikov 1995). Both lineal' fits for spectral types earlier and later 
than AO are overplotted. 
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APPENDIX A: SPECTRAL TYPE - MASS RELATION 

In order to determine the masses of mai n sequen c e star s, we 
used the Stellar Mass Catalogue (SMaC) of iBelikovl i ll 9951) that 
contains the dynamical determination of masses of binary systems 
of luminosity class V, as well as their spectral types. In Fie. IA H are 
plotted the masses of the stars of this catalogue as a function of 
their spectral type. We observe that two linear trends can be drawn 
below and above spectral type AO. We performed a linear fit to both 
trends and obtained the following formulae : 

- for spectral types earlier than AO 
/ M \ 

log — = 1.9 - 0.078 X spt (Al) 



for spectral types later than AO 
I M\ 

log — U 0.70 - 0.016 x spt 



(A2) 



where spt increases by 1 for every spectral subtype, beginning at 
for spectral type O0. Thus for spectral type B0 spt is 10, for AO it 
is 20, and so on. 



APPENDIX B: MASS - RADIUS RELATION 

In order to determine the radius of main sequence stars, we 
used the Catalogue of orbital elements, mass and lu minosities of 
close double stars o flSvechnikov & Bessonova] dl984h that contains 
the mass and radius determinations of both components of eclipsing 
binary systems of luminosity class V. In Fig. lBll are plotted the radii 
of the stars of this catalogue as a function of their masses for masses 
larger than 1 M G . We observe a linear trend at the lower limit of 
the data set, and a scattering that increases towards low mass. We 
suggest that the data points that are situated well above the general 
trend, have a wrong published luminosity class. These points are 
considered as outsiders, and have not been included into the fit. We 
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Figure Bl. Radii as a function of masses for the stars of the catalogue of 
I Svechnikov & Bessonova 1984). The linear fit is overplotted. The rejected 
data points are plotted in light grey. 

performed a linear fit to the remaining data points and obtained the 
following formula : 

log|Aj = 0.021 +0.074* log |-^J (Bl) 



